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Abstract. We consider a system formed by an array of Bose-Einstein condensates 
trapped in a harmonic potential with a superimposed periodic optical potential. 
Starting from the boson field Hamiltonian, appropriate to describe dilute gas of bosonic 
atoms, we reformulate the system dynamics within the Bose-Hubbard model picture. 
Then we analyse the effective dynamics of the system when the optical potential depth 
is suddenly varied according to a procedure applied in many of the recent experiments 
on superffuid-Mott transition in Bose-Einstein condensates. 

Initially the condensates' array generated in a weak optical potential is assumed to 
be in the superfluid ground-state which is well described in terms of coherent states. 
At a given time, the optical potential depth is suddenly increased and, after a waiting 
time, it is quickly decreased so that the initial depth is restored. We compute the 
system-state evolution and show that the potential jump brings on an excitation of the 
system, incorporated in the final condensate wave functions, whose effects are analysed 
in terms of two-site correlation functions and of on-site population oscillations. Also 
we show how a too long waiting time can destroy completely the coherence of the final 
state making it unobservable. 
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1. Introduction 

Bose-Einstein condensation, originally observed in a dilute atomic gas trapped in a 
harmonic potential PjJ |2] , is today obtained in a variety of experimental configurations. 
Experimental efforts have allowed to realize setup in which condensates are achieved into 
one-, two- or three-dimensional optical lattices [SHHElEl, that is arrays of microscopic 
potentials induced by ac Stark effect due to interfering laser beams. Such microscopic 
potentials are often superimposed to the trapping harmonic one and give rise to a 
fragmentation of the condensate. 

In the very recent experiments on Bose-Einstein condensates (BECs) in optical 
lattices, dynamically active states have been generated either by accelerating or by tilting 
the optical lattice or by shifting the harmonic potential trap [3J 0J OH E] Nevertheless, 
in such experimental realizations essentially classical/superfluid regimes have been 
explored and the corresponding dynamics results to be quite well described by the 
discrete Gross-Pitaevskii equation (GPE). Phenomena as Bloch oscillation, nonlinear 
Landau-Zener tunneling, Josephson junction current, can be explained in terms of the 
band structure entailed by the GPE in the spirit of solid-state physics jZj . 

The opposite regime, where a low number of bosons per well or a strong optical 
potential require a quantum description of the dynamics, has been recently explored in 
some experiments [HI E] • A strong interplay between quantum and classical regime takes 
place in such experiments when, for example, the quantum phase transition from the 
superfluid to the Mott insulator regime is generated, or when the collapse and revival 
of the bosonic wave functions is observed. 

In the present paper we consider an experimentally realistic system constituted by 
a dilute gas of N ultracold bosonic atoms, trapped in a harmonic potential and loaded 
into an one-dimensional optical lattice of M wells. At the beginning, we derive an 
effective dynamics by reformulating the second-quantized many-body Hamiltonian, that 
well describes the dynamics of this system, within a generalised Bose-Hubbard model 
(BHM) picture. In such a way we know the effective Hamiltonian dynamical parameters 
as a function of the microscopic system constants as well as of external variables such 
as the magnetic-optical potential strength. Based on such effective picture, we study 
the system dynamics when the optical potential depth is quickly varied. 

Initially, we consider the situation where the system is in the ground-state involved 
by a weak optical potential. Clearly, the discrete GPE should be applied in this regime to 
recognize the ground-state configuration jTUl HH E2| • At time t = the lattice potential 
depth is suddenly increased, so that the tunneling amplitude between neighbouring 
wells quickly drops to zero. Consequently, the system enters in the Mott regime in 
which the time evolution requires a quantum description. After a waiting time t', the 
optical potential depth is suddenly decreased to the initial value. From this time on, 
the discrete GPE gives a satisfactory description of the system time evolution. The 
initial conditions for the following mean-field evolution driven by the GPE thus stem 
from the quantum state emerging from the Mott regime. We compute the quantum 
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state describing the system as a function of time and of the Hamiltonian parameters, 
and show that the phase shift between the condensates of neighbouring wells exhibits a 
strong time dependence while the well populations undergo oscillations. The system is 
thus taken in an excited state f3J El 113 • Further, we show both that the final state 
coherence dramatically depends on the waiting time t', and that the site wave-function 
phase coherence is destroyed when t! is increased. 



2. Space-mode approximation. 

The Hamiltonian operator for a dilute gas of bosonic atoms in a harmonic trapping 
potential Vff(r) = H 3 =1 mQ 2 r 2 /2 with the additional one-dimensional optical lattice 
potential Vl(v) = h 2 to 2 sm 2 (kri) / (AE r ) , [k is the laser mode and E r = h 2 k 2 /(2m) is the 
recoil energy) has the following form 

H = Jd 3 rt[) + (r) 

where ip(r) (tp + (r)) is the annihilation (creation) boson- field operator for the atoms 
in a given internal state, a s is the s-wave scattering length and m is the atomic mass. 
The space- mode approximation jlfij . which allows us to reformulate the system dynamics 
within the BHM picture, is performed as follows. Let Vj be the parabolic approximation 
to V = Vl + T, 3 =2 mQ 2 r 2 /2 in = (jn /k,0,0) the locations of V local minima. We 
assume the energies involved in the system dynamics to be small compared to the 
excitations of the single well ground-state. Thus we can expand the boson field operators 
in terms of Wannier functions ip(r,t) = J2jU*(r)dj(t). In the last equation j runs on 
the optical lattice sites, Uj is the single-particle ground-state mode of Vj with energy 
eigenvalue e(u) = h{oj + Q 2 + ^3)/2. By substituting the previous expression of ip(r,t) 
in the previous Hamiltonian and keeping the lowest order in the overlap between the 
single-well modes, we find the Bose-Hubbard Hamiltonian 

T 

H = J2[Urii(n t - 1) + Xi-rn] - (aj aj + h.c.) , (1) 

i <ij> 

where the operators rii = afai count the number of bosons at the i-site of the 
lattice and the annihilation and creation operators aj and satisfy the standard 
commutation relations [a», a+] = 5ij. In Hamiltonian ^ parameters are defined as 
follows. U := a s Qo^JmhQ / (2n) is the strength of the on-site repulsion, in which we 
have set f2 = tftuTl^h. The site external potential is Xj := e(u) + j 2 Tc 2 h 2 {lf/(4E r ), 
and T := —2 J d 3 rUj[V — Vj±i]uj±i is the tunneling amplitude between neighbouring 
sites. The indices i,j (E Z label the local minima Xj = irj/k of V throughout the lattice 
and Vj = muj 2 (ri — Xj) 2 /2 + Y? l=2 mVL 2 r 2 / 2 . The total number of bosons iV = T^jUj is 
conserved. ^From now on we shall consider a gas of repulsive atoms, thus U > 0. 

BHM ((H) (T21 HH1 IIHI 120] was introduced as model for superconducting films, 
granular superconductors, short-length superconductors, and arrays of Josephson 



-— V 2 + Vff(r) + V L (r) 
2m 



4irh a s 
2m 



d 3 ri[j + (r)^ + (r)^ (r)-?/>(r) , 
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junctions [2H 122] • More recently, some authors (2HI I2H I2H] suggested to describe the 
dynamics of ultracold dilute gas of bosonic atoms trapped in an optical lattices by 
means of a BHM. Experimental results have shown that the essential physics of arrays 
of coupled BECs is captured by BHM (HI El- At zero temperature, the ground state of the 
homogeneous version (Xj = const) of system described by Hamiltonian ((TJ) undergoes 
a quantum phase transition from the superfluid (SF) phase to the Mott insulator (MI) 
one [13 HHJ EH 1201 121 1221- For values of T/U strong enough the ground-state of © 
is a SF and it is well described by a wave function exhibiting a site independent phase 
CH CS] . When the lattice potential depth u is increased and, correspondingly, T/U 
is decreased the system ground-state can manifest two behaviours. If the total number 
of atoms N is commensurate with the site number M the system ground-state is a MI 
with vanishing global compressibility, otherwise it is a SF. In the case of inhomogeneous 
potentials, as that resulting from the confining trap in system (JIJ, there exist Mott 
insulating regions for T/U below a threshold even without commensurate filling [TU| l26j. 

3. Superfluid initial state. 

We consider a set of initial conditions in which the Hamiltonian parameters entail a 
superfluid ground-state. These are achieved in the limit where T/(NU) » 1. For 
N large and T/(NU) over a suitable threshold, it is widely accepted that the low- 
temperature dynamics of can be described by a discrete version of GPE 
This semiclassical limit can be accomplished recalling that when the tunneling term 
dominates the on-site repulsion one (namely T/(UN) » 1), the Glauber coherent 
states give rise to effective solutions for the quantum problem entailed by Hamiltonian 

H3 HOI HE] • Thus, a reasonable solution for the BHM ground-state in the regime of 
interest, can be obtained within a coherent state variational picture based on applying 
a time-dependent variational principle on coherent-state trial state. By means of this 
procedure (for details see [2311011201), the quantum dynamics generated by Hamiltonian 

can be reformulated in terms of a classical dynamics generated by an effective 
Hamiltonian 7i. Hence, following this procedure, we assume the system dynamics to 
be described by the trial state |\&) = exp(iS/h)\Z) , where \Z) := Hi\zi) (see [101 1201) i s 
written in terms of Glauber coherent states as 



(recall that a, |0) = 0, and that their defining equation is a.i\Zi) = Zi\z/) with 
Zi G C). The effective equations of motion are achieved by a variational principle 
from the effective action S = J dt[iHj(zjZ* — z*jZj)/2 — 7i], associated to the classical 
Hamiltonian T~t(Z, Z*) := (Z\H\Z), through a variation respect to Zj and z\. Hence the 
time-dependent trial-state parameters Zj = (^f\aj\^) represent the classical canonical 
variables of the effective Hamiltonian dynamics and satisfy to the Poisson brackets 



Zi) :=e-^' 2 £^(a!) n |0> 
Z^n n! 



n=0 



(2) 
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{z*,Z(} = i5jg/h. After some algebra we obtain the classical Hamiltonian 

n = Y,[u\z j \* + \ j \z j \ 2 -Z( z * Zj+1 + c .c. )], (3) 

j 

where j runs on the chain sites: j G I M with I M = {0, ±1, . . . , (M — l)/2} or 
J M = {±1/2, ±3/2, . . . , M/2} when M is odd or even, respectively. The related 
equations of motion are the following 

ihzj = (2U\zj\ 2 + Xj)zj - \{zj-\ + Zj+i) , (4) 

with j G I M , and with the complex conjugate equations. 

The ground state of Hamiltonian (jHJ) is determined by studying the variation of 
7i — xTV, where the Lagrange multiplier x has been introduced to explicitly incorporate 
the conserved quantity M = Y,j\zj\ 2 . In the limit where the on-site chemical potential 
Xj is slowly varying with the site index j, namely \Xj — Xj + i\/U << 1, we have 

+ Zj + i « 2zj. An approximate solution for the SF configuration is thus given 

by 



M' X = 2UN + M'X-M'T, Zj = ^ - Qi^A e <* , (5) 

where A = S |j|g/ Xj/M' and M' = min(M,q) is determined by finding out the 
maximum integer q such that 2UN ± q(X — X q ) > 0. This solution represents the 
discrete version for the Thomas- Fermi approximation [2H] • The corresponding energy is 

E gs = N/a[l - M't - a(X + o/A{5X) 2 )} , (6) 

in which we have set M'(5X) 2 := S j£/ / (Xj - A) 2 , a := M'/(UN) and r := T/(UN). 
Concerning the three coupled condensates system within the SF regime, a thorough 
study has been made in Ref. [2H] . 



4. Superfluid to Mott-insulator transition. 

At the time t — the optical-lattice potential is suddenly increased. This is achieved 
by varying the potential intensity according to, for example, a tilted slope: uj(t) = 
ou[l+(w —l)t / Tf,], where r& is the time scale for the jump and w is the amplification factor. 
Consequently, the tunneling amplitude goes to zero T[o;(t)] — ► as an exponential. In 
fact, standing the definition T := —2 J d 3 ruj[V — Vj±x)uj±i, where Uj is the harmonic 
oscillator ground state involved by the quadratic single particle potential Vj, by direct 
analytical calculations we get 



Also U and Xj are modified when changing the optical potential amplitude, but their 
dependence on u(t) is much less dramatic. In fact we have U[uj(t)] = U[uj(t) / 'uj] 1 / 2 and 
Aj[a>(t)] = Xj + K[u(t) — u)]/2. While the potential depth is changed, < t < u, we 
suppose the state of the system to be not modified. 



T[u{t)\ 



tiW(t) 

4E r 



ir 2 2E r 

1 H — 

2 hu(t) 
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To apply the sudden approximation, the time scale r b characterising the potential- 
depth jump must be fast compared with the tunneling time between neighbouring 
wells, but slow enough to prevent the condensate excitations in each well, namely 
2h/oj « r b « h/T (where T is the larger hopping amplitude). In fact, doing so 
the system persists in the lowest band and the effects due to the hopping term result 
negligible. Thus, a general initial state D{n} c n| n ) ; will have the straightforward to 
compute time evolution J2{ n } exp{— z/fr J2j[Un^ + Xjrij]}c a \n) . Concerning the sudden 
approximation, we recall that it is usually used for calculating transition probabilities in 
the case when the Hamiltonian changes rapidly within a short time interval (presently 
this is identified by t = and t — T b > 0). One simply assumes the reaction of the 
initial state to the quick Hamiltonian change to be negligible. So one can approximate 
the transition amplitude by assuming: (out\U(r b ,0)\in) ps (out\in). Such an argument 
works only if impulsive forces are absent, which, otherwise, could generate finite states 
change even if applied for infinitesimally long time. In our system such kind of forces 
are absent. In fact, in the time interval (0, r&), where the Hamiltonian parameters are 
time dependent and the jump of the optical-potential depth is driven by w, by means 
of the coherent state representation of the path integral, we have (SOj 

i r b 

h 



(z\U(r b ,0)\z) = I V[z]exp 

where 



dtL(t) 



(8) 



L (t) = |£ - m^{t)]-n[z{t) : z*{t)^ 

is the Lagrangian of the effective path- integral action. Since L{t) can be shown to have 
no singular behaviours as r b — > 0, namely lim n ^ + J Tfc dtL(t) = 0, the above formula 
implies that the dynamical evolution of the initial state is driven by r&. So the shorter 
Tf, implies the smaller changes of U(r b , 0)\z). 

In the Mott regime {T/U << 1), namely after the potential jump (t > r b , we 
will assume r b = 0), the classical description of dynamics is no longer valid and the 
Schrodinger equation is necessary to describe the time evolution. The latter is generated 
by Hamiltonian (0) where T = 0, U = U{wuo) =: U and Xj = Xj(wu) =: Xj are assumed. 
As assumed above, the system initial state is described by the coherent state (J2J) whose 
the quantum evolution is given by 

|(f)> = e-* ffi II l*> =n^E ^ V f§ ni e*S"®\TH) , (9) 



ieI M i£/ M ni=0 

where we have set Vi{t) := exp[i/h(U + Xi)t] and u(t) = Ut/%. The quantum time 
evolution does not preserve the coherent state structure of state (J2J), in fact, the term 
exp[— infu(t)] in © breaks the coherent state form. By direct calculations, it is easy to 
show that the quantum evolution in the Mott-regime time interval entails the following 
expectation values 

((f)K + Gil(f)) = N 2 , 
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Zj (t) := ({t)\aj\(t)) = ^•exp||-A j t-z|z i | 2 sin 2 [u(t)]|exp{-2|^| 2 sin 2 [u(t)]} , (10) 
((t)\a} +iaj \(t)) = z* +1 ^exp{-2(| % | 2 + \z, +1 \ 2 ) sm 2 [u(t)}} x 

xexp{~(\- X j+ i)t - i{\zj\ 2 - \z j+ i\ 2 )sm 2 [u{t)}} . (11) 

Such equations display that, during the quantum time evolution, the wells population 
does not change, whereas the site wave-functions Zj(t) are dynamically active and driven 
by more then one characteristic times. The modulus of Zj(t) is a periodic function of t 
\zj(t)\ = \zj\ exp{— 2\zj\ 2 sm 2 (Ut/h)} with period T m = nh/U. The phase of the site 



0.2 0.4 0.6 0.8 

time 



1.2 1.4 



Figure 1. This figure shows the time-evolution of the modulus of the central-site 
wave- function Eo(t) = z o(t)/ z Q as given by Eq. i|l(J|) . It displays a periodic behavior 
whose period is T m — ttTi/U w 0.396 sec. Notice the "almost impulsive" periodic 
behavior of £o(t) which is mostly almost vanishing. 



wave-functions tpj := arg[^-(t)/|^-(t)|] = Xjt/h — \zj\ 2 sin 2 [u(t)] is driven by T m , the 
period of sin 2 (/u(t)], and by the = 2nh/Xj the site dependent periods involved by 
the external potential Xj. Furthermore, in each site j where \zj\ 2 exceeds the value 2ti 
the condition \zj\ 2 sm 2 [u(t + T z )] = \zj\ 2 sin 2 [u(t)} + 2tt implies a further characteristic 
time T Z] for j e hi- 

Also, the site-dependent external potentials Xj induce a dephasing between the 
Zj(t) and zj + i(t) representing the condensate states (namely the site wave function) at 
sites j and j + 1. In fact, it results ipj — fj+i = —(nhQ 1 /2) 2 (2j + l)t/{TiE r ) — i\zj\ 2 — 
\zj + i\ 2 ) sin 2 [w(t)] that shows as the dephasing increases along the lattice. The difference 
between the phases of the condensates announces that the system is no more in the 
ground-state O EH Hi] and give rise to excited configurations. 

Figures ^01 have been achieved by considering e realistic experimental configuration 
with 10 4 85 Rb atoms, a harmonic trapping potential frequency Qj pa 50 Hz (j=l,2,3), 
and a laser mode k ~ 10 7 m. The optical potential amplitude, initially of the order E r , 
is suddenly increased to 30E r . 
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Figure 2. Trajectories of the site wave-functions £j(£) = Zj(t)/zj in the complex 
plane (abscissae and ordinates refer to real and imaginary part, respectively) for 
< t < AT m . Left and right figures refer to site j = (the central one) and site 
j = M/3, respectively. Notice that, unlike T, (t), T, M / 3 (t) is not periodic. 
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Figure 3. The evolution of the spatial correlation of adjacent sites is represented by 
the parametric plot of = ((t)|a~ f +1 a J ;\(t))/z* +1 Zj, with j = M/3 and < t < 4T m . 
A strong dephasing between the sites under concern is evident. 



5. Mott-insulator to superfluid transition. 

Times t > t' correspond to the third stage of the dynamics. The optical-potential 
depth is quickly decreased (in a time of order u w 0) to the original value u thus 
restoring the regime with T / {NU) » 1. In this case an approximate description of the 
system dynamics within the semiclassical variational picture applied in the superfluid 
regime should be again applicable. The initial conditions for the third stage of the 
system evolution can be easily shown to be represented by a superposition of Glauber's 
coherent states. In fact, the initial state is obtained from Eq. (jHJ) by setting t = t' 

id m nj=° \J n r 

Here £j = exp{ — \zj | 2 /2}, z'j = ZjVj(t'), and u' = u(t'). By using the identity 
lim / dxexp[— (p + e)x 2 — inx] = exp[— n 2 / (Ap)]J it / p , 
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with p = —ij (4m') , this state can be written in a very suggestive form as a superposition 
of coherent states. Direct calculations give 

ko)= n r ^e-^v^\>e-^) = n r^^uoi^), (12) 

■ j J-00 Ayixu J-00 

where the state labeled by z'J = Zje~ lXj is the normalized coherent state given in Eq. 
(J2J) with Zj = z". If one expresses \(t')) as \(t')) = IT/ei the new trial state 

accounting for the evolution after the second potential-depth change, might be expressed 
as 

!(£;)>= II ^)l(Oi> 

with D(£j) = exp(^ja+ — where the time behaviour of new dynamical parameters 

£j, £* occurring in the exponential terms (actually these are coherent-state displacement 
operators) must be reconstructed by implementing once more the time-dependent 
variational procedure. Due to the properties characterising the the displacement- 
operators action on a coherent state \z) (D(£)\z) = exp[ilm(z* (,)}\z + £)) the final 
form of the trial state is 

/oo 

where the kernel K(xj,u') is defined implicitly, and £j(t) = at t = and <pj = 



6. Final discussion. 



In the present paper we have considered a performable experimental process exhibiting 
a strong interplay between classical (SF) and quantum (MI) regimes. We have described 
the dynamics of an array of BECs when the optical potential depth is quickly varied. The 
process we have considered forces the system to go through an intermediate quantum 
regime. As a consequence of this, the system loses its semiclassical character assuming 
the form of an excited state that cannot be represented as a simple direct product 
of coherent states. Eqs. (jllj) show that collapsing/revival phenomena occur whose 
characteristic time scales have been recognized. Further, the presence of the harmonic 
external potential appears to responsible for a strong site dephasing. When the optical 
potential depth is lowered again, the resulting state has been shown, within the previous 
section, to be a superposition of coherent states represented by the integrals of eq. (fT2")) . 
In performing the integration on Xj, at each site, each coherent state \zje~ lXj ) contributes 
with a phase e _iXj . The latter might have a destructive effect for increasing vl when 
calculating the expectation value of the physically relevant operators of the model. This 
can be seen by re-expressing the integrals in (|T2"|) in the form 

dx ;„1IU».I\ „„„r^ .--.w„tl 



'iru' 



D ix /(Au 1 ) exp[(— ix)z'a<] 
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Also, since the term exp[ix 2 /(4m')]/ y/u' is rapidly oscillating outside the interval 
[—y/nu', y/nu'], the major contributions is expected to come from the integration 
on this interval. Observing that the corresponding coherent-states phases change in 
— y/ixu' < x < y/nu', the simplest way to reduce the decoherence effects may be achieved 
by imposing vl << 1, that is t' « %/U. In view of their complexity this problem and 
the evaluation of the time behaviour of the state emerging from the Mott 

regime will be discussed in a separate paper. 
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